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Abstract 

Markov chain approximations of symmetric jump processes are investigated. Tight- 
ness results and a central limit theorem are established. Moreover, given the gen- 
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1 Introduction 



Let Y = (Y t ) t be a time-continuous Markov chain on Z d . It is a natural question whether 
the sequence (Y n ) of Markov chains defined by Y t n = rr x Y n0lt , a G (0,2], tends to some 
reasonable process for n — > oo. The case a = 2 is known as diffusive scaling and leads to a 
diffusion process under certain assumptions on (Y t )t, see the classical Donsker's Invariance 
Principle of |Don51] for the Brownian motion and chapter 11 of [SV79] for diffusion pro- 
cesses in non-divergence form. In the case of symmetric processes Stroock and Zheng derive 
in [SZ97] a central limit theorem for continuous-time symmetric Markov chains of bounded 
range. In a recent paper, Bass and Kumagai [BK06J remove the restriction of bounded 

range by replacing it by a second moment condition. In both publications, the generator 

d 

of the limit object is of the form Lu(x) = Yl 9 Xi (aij(-)d Xj ) and a formula is provided 

how the diffusion coefficient functions djj(-) can be computed from the conductivities of 
the chain (Y t ). The other direction, i.e. constructing a sequence of approximating Markov 
chains for a given diffusion matrix is not less important and one of the main results of 
[SZ97j . 

The aim of this work is to prove results analogous to ones of |SZ97j . [BK0 6] in the case 
where the limit object is a symmetric jump process with corresponding Dirichlet form 
(£,D{£)) given by 

£(f,g) = l2 J J {f(y)-f(x))(g( y y)-g(x))k( y x,y)dxdy, k(x,y) = k(y,x) , 

(1-1) (R d xR d )\diag 

D{£) = Cjm £ \ where £ 1 {f, f) = £{f, f) + \\f\\ I?1 
and generator L given by 

^,-5. / (*>-«,,>*.„*. 

\y-x\>e 

Therefore, we study Markov chain approximations for a certain class of symmetric jump 
processes. In [SZ97] . [ BK0 6] the generator of the limit object is a uniformly elliptic op- 
erator. In our situation the equivalent concept of uniform ellipticity would be given by 
k(x, y) > c\x — y\~ d ~ a V \x — y\ < r for some c > 0, r > 0, a G (0, 2). One feature of our 
approach is that our central limit theorem allows for cases where such an estimate does 
not hold, i.e. the limit process may be a pure jump process which is anisotropic in some 
sense. The level of anisotropy is limited since our approach uses a-priori bounds for the 
modulus of continuity of the heat kernel. As discussed in [B BCK06] these bounds fail for 
very irregular jump measures. 
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There are several other contributions to the question how to approximate Hunt processes 
given by Dirichlet forms, see |MRZ98j . [MRSOO] and the references therein. However, our 
results are not covered by these works. We close the introduction by commenting on the 
differences between this work and [BK06J, [SZ97J. 

• The limit object in [SZ97J and |BK 06] is a diffusion whereas here it is a jump process. 

• The main result of |BK06] is a central limit theorem. In addition to such a result we 
establish an approximation result for a given symmetric jump process. Theorem 12.31 
should be contrasted with Theorem 3.9 from [SZ97J . 

• Our assumption (A5) differs from (A5) of [BK06j significantly. First, we do not 
assume continuity of the coefficients of the limit process. Second, we assume only L] oc - 
convergence of conductivities which is substantially less than uniform convergences 
on compacts. 

The paper is organized as follows. In section [2] we present our assumptions and results. 
We provide a detailed discussion of the assumptions and some definitions and notation. 
Furthermore, an auxiliary result on equivalent norms on the Sobolev space H a ^ 2 (R d ) is 
proved. Sections [3] and H] provide the proofs of Theorem 13.61 and Theorem 14.61 both of 
which are crucial to the proof of our main results. In section we prove Theorems 12.11 and 
12.21 Theorem 12.31 is proved in section El 

Acknowledgments: The authors express their gratitude to Professor Z.-Q. Chen and to 
Professor T. Kumagai for helpful comments on the presentation of the results. 

2 Assumptions and results 

We formulate our assumptions and results in section 12.11 In section 12.21 we provide a 
detailed discussion of the assumptions. Section 12.31 is devoted to a result on equivalent 
norms on H a / 2 (R d ), a e (0,2). In section [231 we define and list various further objects 
that we deal with in this article. 

We denote the counting measure by p and the Lebesgue measure on R d by A. In our context 
we also deal with the function spaces L 2 \p~ lr L d , p~ d p) where p > 0. The scaling factor p~ d 
in front of p is natural from a geometric point of view. Write B p (x,r) := B(x,r) D p _1 Z d 
for the r-ball around x in p~ lr L d . We also use the notation p p = p~ d p. For x G R d we use 
the abbreviation I |. For iGRwe write \x~] instead of max{Z 6 Z : / < x}. 
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For a point x G R d we denote by [x] n the element of n Xr L d satisfying ([x] n )i = n x \riXj\ 



for alii = 1, . . . , d. 



2.1 Formulation of assumptions and results 

Let (C n ) ne iN be a sequence of conductivity functions C n : n~ lr L d x n~ lr L d — > [0,oo). The 
following assumptions will be important. 

(Al) C n (x, y) = C n (y, x) and C n (x, x) = for all n G IN, x, y G ri _1 Z d . 
(A2) There exists k\ > such that for all n G IN, x ^ y 

C n (x,y) <Ki\x- y\~ ~ a . 

(A3) There exist Nq G IN and K2 > with the following property: For any n G M, 
i,l/6 n~ 1 Z d , x y there are elements Zq\ . . . , zf G n _1 Z d , Z < jV , Zq = x ' 
= 1/ satisfying 

C^,^)^^^-^, i = 0,...,Z-l. 

Moreover, for any pair (CO G n~ lr E d x n _1 Z the number of pairs (x, y) G ri _1 Z x 
n _1 Z d such that £ = z^'^ and £ = ^fc^ for a k is bounded by A" . 

For given x, y G n _1 Z d , i^j/ we call the ordered set {zq 2 ^, . . . , 2j ^} above a chain and 
/ the length of the chain. The above assumptions are essential for our approach. For a 
mere technical reason discussed below in detail we need an additional assumption: 

(A4) There exist 0i > and k 3 > such that for all x G n~ 1 1t and r > ®\n~ l 
fx({y G B n (x,r) : C n (x,y) > k 3 \x - y]'^ }) > § //(£ n (z, r)) . 

It is important for our results that the constants Ki, K2, K3, Nq, Oi appearing in (Al) 
through (A4) do not depend on n G IN. We associate to C n a symmetric discrete-time 
Markov chain X n = (X£) fceM by 

(2.1) P-W = »)= ^b) - 

Let F n = (V < n )i be the symmetric time-continuous Markov chain that has the same jumps 
as (Xh) while its holding time in the point x is exponentially distributed with parameter 
J2zen- 1 z d C n (x, z). Note that each Y n starting in x G n -1 Z d corresponds to a probability 
measure on D([0, 00); R d ), the space of right-continuous paths in R d having left limits, see 
[EK86J, [Bil99] for properties of D([0, 00); H d ). Our first result reads as follows: 
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Theorem 2.1 Let {C n ) n be a sequence of conductivity functions satisfying (Al) through 
(A4)- For x n G n~ lr L d , x n — > x G R d the laws ofY n starting in x n are tight in D([0, t ]; H d ) 
for any t > 0. 

For a reformulation and the proof of this result see Theorem 15.11 

In order to establish a central limit theorem one needs to prescribe the behavior of C n 
for n tending to infinity. For x G n~ lr L d set Q. n {x) = Ylf = x[xi,Xi + and £2 n = 

U { a. ( •<•)}. 

(A5) There exists a measurable function k: R d x R d — > [0, oo) such that for any compact 
subset K of R d x R d \ { (x, x) : x G R d } the funct ions (x, y) h- > C"([x] n , [y] n ) converge 
in L}(K) to •) for n — > oo. 



(B) There exist M G IN and A 2 > with the following property: For any e > 0, n G IN 
and O, Q G n there are elements P<f ' S) , . . . , P; ( °' S) G0 n ,K M with P^' S) = O, 
Vi°' Q ^ = Q and satisfying 

k(x,y)t{\ x - y \> £ }dxdy > A 2 // l{| x _ y |> £ }|a; - y|"^ Q rfxrfy Vj < I - 1 . 



'j X/ 3+1 

Moreover, for any n G IN and any pair (1Z,S) G 0. n x n the number of pairs 
(0, Q) G Q„ x £} n such that ^ = vf' Q) and 5 = ^J S) for a A; is bounded by M . 

Again, we call the ordered set {Vq°' Q \ . . . above a chain and / the length of the 

chain. Although, to some extent, (B) is a continuous analog of (A3) it does not follow 
from (A3) and (A5). Such an implication could easily be achieved by adding an additional 
assumption. In order not to weaken Theorem 12. II we prefer to work with (B) separately. 

Here is our central limit theorem. 

Theorem 2.2 Let (C n ) n be a sequence of conductivity functions satisfying (Al) through 
(A 5) and (B). Let 3€ be the symmetric jump process associated to the regular Dirichlet 
form (£,D(£)) given in ( SUP and Af the properly exceptional set. Then, for x n G n~ lr L d , 
x n — » x G R d \ Af the laws ofY n starting in x n converge weakly in D([0, i ]; R d ) to the law 
of S£ starting in x. 

Remark: The assumptions (A3) and (B) are technically involved and cover anisotropic 
situations. In fact, (A3) and (B) are trivially satisfied in the isotropic case, i.e. if C n 
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satisfies (Al), (A2), (A5) and C n (x,y) > c\x — y\ d 01 for all n G IN and \x — y\ > n l K 
for some K > 0, c > 0. Even in this case our theorem is still interesting and new. 

It is necessary to allow for some exceptional set in Theorem 12.21 However, due to results 
in [CK03J the set M can be removed or assumed to be empty in several situations. There 
is no need for an exceptional set in our third result, Theorem 12.31 As in Theorem 3.14 of 
[SZ97j we give an explicit construction of approximating Markov chains. 

Theorem 2.3 Let k: R d x R d — > (0, oo) be a measurable function satisfying k(x,y) = 
k(y,x) and 

(2.2) k 4 \x — y\~ d ~ a < k(x, y) < k 5 \x — y\~ d ~ a 

for all x,y G R d ; x ^ y with some positive constants < k&. Define the conductivity 
functions C n : n~ lr L d x n~ lr L d — > [0, oo) by 

{0 for \x - yloo < n' 1 , 

n 2d f KZ,0dZd( for\x-y\ OQ >2n- 1 
|x-CU<n-V2 
|2/-Cloo<™-72 

Let 3£ be the Hunt process corresponding to the Dirichlet form (£,D(S)) given by /{!.!]) . 
Then the sequence of processes corresponding to C n converges in the sense of Theorem \2.2\ 
to 3£ for any starting point x G R d . 

Remark: As becomes clear by the discussion below (A4) allows for quite general cases of 
sequences C n . In addition to it, in light of Lemma I2~4l and Lemma [3TB1 it is very likely that 
(A4) can be dropped. This would imply the possibility to weaken the lower bound on k 
assumed in (12.21) substantially. 

2.2 Discussion of assumptions 

We illustrate assumptions (Al) through (A5) introduced in Section |2~T| First, let us look 
at (Al) through (A4). If, for a fixed scale ne^,C n : n~ Xr L d x n~ lr L d -> R + satisfies (Al) 
through (A4) then the same holds for the conductivity function C : Z d x Z d — > R + with 
C(x, y) = n^^^CirT^x^ n~ x y) for x, y G Z d with the same constants d, k±, K2, N where the 
chains in (A3) have to be scaled in an obvious way. Naturally, the interaction radius Oi in 
(A4) is multiplied by n. Since C is the appropriate conductivity function corresponding to 
the process Y n scaled on Z d in the obvious (" a-stable" ) way, it is sufficient to understand 
(Al) through (A4) for a single conductivity function C on Z d x Z d . In addition, the 
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main results of Section 3 and 4 are scale-invariant, i.e. the constants appearing are scale- 
invariant, and depend only on the constants in the assumptions. Therefore it again suffices 
to prove them for a fixed conductivity function. 

(A3) and (A4) are stable under perturbations of the conductivity function near the diag- 
onal. This reflects the fact that in our central limit theorem jumps smaller than a fixed 
R > have no influence on the limit process. Let C and C be conductivity functions 
such that there exists R > with C(x,y) = C(x,y) for \x — y\ > R. If C satisfies one 
of the assumptions (A3), (A4) then the same assumption also holds for C(x,y) with the 
same constants. On the other hand, if C n ,C n : n~ lr L d x n~ 1 Z d — > R + are to sequences 
of conductivities with C n (x,y) = C n (x,y) whenever \x — y\ > Rn^ 1 and if (C n ) satisfies 
(A5), then (A5) also holds for (C n ) with the same limit function k. 

(Al) implies the process to be symmetric while (A2) bounds C(x,y) from above by the 
conductivities of a rotationally symmetric a-stable Markov chain on Z d . (A2) gives in 
particular 

sup ^ C ( x iV) < °°- 

(A3) is much more technical. It implies a certain kind of irreducibility of the associated 
Markov chain. Additionally, it takes into account the highly non-local nature of our objects. 
Roughly it says that every two points x, y can be connected by chaining together at the 
utmost N Q jumps where the probability of each jump is bounded from below by a constant 
multiple of \x — y\~ d ~ a while at the same time one has enough of these connecting jumps. 

(A2) and (A3) imply together \z\ x ' v) - z%f\ < («i/ K 2 ) 1/id+a) \x - y\. This leads us to the 
following necessary condition for (A2) and (A3). 

Lemma 2.4 Assume (A2) and (A3). Then there exist 7 e (0,1), ©i > and k 3 > 
depending only on Ki, k 2 , N , d and a such that for all x e , r > ©i 

(2.3) ^{{y G B 1 (x, r) : C(x, y) >k 3 \x- y^}) > l^B^x, r)) . 

In particular, if the conductivities are stationary, i.e. C(x,y) = C(x — y) then 
fx({h E B\0,r) : C(h) > ^1^}) > 7M^(0,r)). 

Proof: First notice that (A2) and (A3) imply the existence of c\ = ci(«i, k, 2 , d, a, Nq) > 1 
such that for all I < N , f , C G % d 

l^-el^ciie-ci- 
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Assume r large enough and x £ Z d . Then M = {z^ £ Z d : y £ B 1 (x,r/ci)^ C i? 1 (a;, r). 
By the second part of (A3) 

(2.4) MM) > l^B^r/c^/No] > C j^-^(B\x,r)) . 



Next, set 



Trivially, /i(M) > C3fx[B 1 (x, r)j where C3 = 4e d^y o depends on all constants that appeared so 

far. Assume r > ( 2 c^n ) * ^ an< ^ z ^ Then there is y £ -B(x, r/ci) with z = z[ x,v ^ £ M 
and 

C(x, 4* > k 2 \x - yV d - a > K 2 4 +a r- d ~ a > K2 (^-f +a)/d \ x - ^p^. 

Setting k 3 = K 2 (-^y d+a ^ d and 7 = C3 the set {y £ B 1 (x,r) : C(x,y) > k 3 \x — y\ d ~ a y 
contains M and satisfies ( 12. 3j) . ■ 



Lemma 12.41 immediately implies that under our assumptions a second moment condition 
as in [ BK06] cannot hold. But (12.31) is not sufficient neither for (A2) nor for (A3). Take 
for instance d — 1 and C(x,y) — \x — if x 7^ y, x — y £ 2Z and C(x,y) = else. 

Then clearly (A3) is not satisfied. 

Let us now provide some examples of conductivity functions satisfying our assumptions. 
If C(x,y) \x — y\ d+a stays bounded between two positive constants then C satisfies (A2), 
(A3) and (A4). Hence all cases of |BL02bj are covered by our conditions. In addition, our 
assumptions allow for cases where there are no jumps in the direction of certain cones. 



Example 2.5 Let V := {(xi,x 2 ) £ Z 2 : \x 2 \ < 7 7 > be a double-cone in Z d . Set 
(2.5) C(x, y) =l v {x- y)g(x, y) \x - y\~ d ~ a 

for x,y £ 1? , x 7^ y where g: Z d x Z d — > [a,b] for some < a < b is a measurable, 
symmetric function. Then these conductivities satisfy (A2) and (A3). If 7 is large enough 
(A4) holds, too. 

In fact, if x = {x\,x 2 ) and y = (yi, y 2 ) set z xy = (xi + [(2 + 7 -1 ) \x — y|] , x 2 ). Then 

x - z *y = (\(2 + 1 - l )\x-y\],0)eV, 
\x-z xy \ < (3 + 7- 1 ) \x-y\ , 
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z x y -y = (x x -y 1 + [(2 + 7 x ) \x - y\] , x 2 - y 2 ), 
\z xy - y\<\x-y\ + V2\(2 + 7- 1 ) \x - y\] < 2(4 + 7- 1 ) \x - y\ . 

Finally, z xy — y G V by 

1 \xi - yi + \(2 + 7 _1 ) \x - J/11 1 > 7f(l + 7 _1 ) \x - y\] > \x - y\ > \x 2 - y 2 \ ■ 

Example 2.6 Define C as in example \2. 5\ with 7 large enough, say 7 > 7/8, and g = 1. 
Set C n (x,y) = n d+a C{nx,ny) ,n G IN", x, y G n~ lr L d . Then {C n ) n satisfies (Al) through 
(A5) with 

(2.6) k(x, y) = l v (x - y) \x - y\~ d ~ a , V := {(xi,x 2 ) G R 2 : \x 2 \ < 7 |xi|} . 



Note also the following counterexample: 



Example 2.7 Set V := {0} x Z U Z x {0} and C(a;, y) = t v (x -y)\x- y\~ d ~ a . Then 
these conductivities do not satisfy (A2) and (A3) since the necessary condition of Lemma 



~2J\ does not hold. 



Finally, let us give the most obvious example of a conductivity function C satisfying (Al) 
through (A5) . 



Example 2.8 Fix a G (0,2). Then the conductivity functions 

C n (x,y) = - -77^—7^ — --lx-2/r , x,yen- lr L d 

v i»> 2 1 ~ a 7r d / 2 T(l - a/2) 1 y| ' y 

satisfy (Al) through (A5). The limit process in the sense of Theorem \2.2\ is the well- 
known rotationally invariant a-stable process. The properly exceptional set M is empty. 



2.3 Equivalent norms on 



So far we have concentrated on a discussion of (Al) through (A5). Let us now look at 
(B). Let k: R d x R d — * [0, 00) be a measurable function satisfying k(x,y) = k(y,x) and 
k(x,y) < Ai \x — y\~ d ~ a for almost all (x,y) with x ^ y and some Ai > 0. In light of 
(Al), (A2) and (A5) this is the structure of kernels appearing in the limit n — > 00. Under 
assumptions (Al) through (A5) there still can be large oscillations of C n in the following 
sense. Fix two sequences of x l n , y l n G n _1 Z d , i G {1, 2}, with \x\— x 2 ^\ — > and \y\-yl\ — > 0. 
Then chains connecting x\ and y\ can be very far apart from chains connecting x 2 n and 
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y 2 , no matter how large n is. Assumption (B) guarantees that this phenomenon can be 
avoided by choosing appropriate chains. Having studied (Al) through (A3) it should be 
clear how to construct examples of kernels k satisfying (B). For instance, the kernel k 
constructed in example 12.61 satisfies (B). 

Let us show that (B) is a natural assumption. Assuming (B) we show that D(£) from (II. ip 
equals H a / 2 (H d ), i.e. (B) determines a class of equivalent norms on H a/2 (R d ), a G (0,2). 
One standard definition of this function space is 

H */2 (Rd) = {fE L 2 {Rd) . mHa/2 1= y/wfw^ + e^fj) < 00} . 

We show how it is possible to replace \x — y\~ d ~ a in the definition of £ a (f,f) by some 
anisotropic kernel k(x,y) without changing the function space. Such a result does not 
seem to be established in the literature on function spaces. 



Theorem 2.9 Let k: H d x H d —>■ [0, 00) be a measurable function satisfying k(x,y) = 
k(y,x) and k(x,y) < Ai \x — y\~ d ~ a for almost all (x,y) with x 7^ y and some Ai > 0, 
a G (0,2). Assume that k satisfies (B). Then there are two positive constants c ,Ci such 
that 

(2.7) c £ a (fJ)<£(fJ)<ci£ a (fJ) V/6«. 

Hence, the regular Dirichlet form (£,D(£)j of II 1.1}) satisfies under (B) 

D{£) = CliWf 1 = H a/2 (R d ) . 



Proof: The second estimate in (12.7ft follows trivially from the upper bound of k. In order 
to establish the first one note 

£*(/,/)= Mm J J (f(y) - f(x))h {lx _ yl > £} \x -y\- d ' a dxdy 

R d xK d 

Denote by zq the center point of a given cube Q G n . For simplicity we assume that the 
length of each chain is equal to M Q . Assumption (B) gets involved in the following way: 

(f(y) ~ f{x)) 2 l{\ x - y \> e} \x - y\~ d ~ a dx dy 



U d xR d 



lim Yl Y (f( z o) - f( z Q)T // 1 {\x-y\>e}\x -y\ d a dxdy 



n— >oo 

oea n se£!„ 



Oxi 

Mo-1 



- M °i™o Y Y Y {f( z r^)-f( z p{°^)) 2 // 1 i\x-y\>e}\x - V\ d °dxdy 
n ^°° OeQ„ Qe£J„ 3=0 " 
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M -l 



-TTrSiSZ E E (^^p ( °; S) ) ~ f(z v (o,Q))) 2 II l {]x _ y] > E} k(x,y)dxdy 
2 n ^°° 0eQn geQ n j=0 



(o,a) vX ,(o,Q) 



< — — lim V V max (f(z (o, e) ) - f{z (.0,0.))) 
Oe£2„ QeQn 

x if \ { \ x _ y \> e] k(x,y)dxdy 

o 
+ 

(M ) 3 



V (0,Q) xV (O.Q) 



< 



A2 " — x 



lim E E (f( z o) - f( z Q)) t {{x _ yl > £} k(x,y)dxdy 



0eO„ Se£!„ 0x 



11 (f(y) - f(x)) Z ^-{\x-y\>e}k(x,y) dxdy. 



A 2 

(R d xR d )\diag 



2.4 Further definitions and notation 

If is a stochastic process and Q a Borel set write t(Q; JT) resp. <r(f2; i£") for the first 
time the process exits resp. enters Q where we omit if there is no danger of confusion. 

Let X = {Xk)k be the symmetric discrete-time Markov chain associated to C by (12.11) . 
A symmetric time-continuous Markov chain Y = (Y t ) t having the same jumps as X can 
be constructed as follows: Take a family (T x ,j) x ew, d jem °f independent random variables, 
independent also of X, such that T x j is exponentially distributed with parameter C x and 
set T xfi = 0. Set Y t = X n for t e E"=o T x 3 ,j , Y%=o T Xj ,j ) • Note that in [BLU2b] the 
holding times of the time-continuous process are exponentially distributed with parameter 
1 leading to different generators and Dirichlet forms. For more details on Markov chains 
we refer the reader to [Nor 9 7] . (A2) and (A3) give uniform bounds on the expected holding 
times of Y. The process Y corresponds to the Dirichlet form 

w*f) = \ E {m - f(y)) 2 C(x,y) . 

x,y£Z d 

with domain D(S) = L 2 (1i d ,fi). We derive properties of Y in the next section by compar- 
ing (S, L 2 (Z d , /i)) to the Dirichlet form (S a , L 2 (Z d , //)) of a rotational invariant a-stable 
process V in Z d , i.e. £ a is defined by 

Uf, /) = ~ E (f^ - f^) 2 \ x - y\~ d ~ a ' « e (0, 2) . 

x,y€l, d 
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By (A2) £(f,f) < ^SMf) for all / G L 2 (Z d ,p). 



Define a family (Y p ) p> o of symmetric time-continuous Markov chains on p lr L d by Yf = 
p~ l Y p a t . Y p corresponds to the Dirichlet form (S p , L 2 (p _1 Z d , p~ d p)) defined by 

(2.8) 

mf, f) = \ E (/(*) - f(y)Y CP ( x > y)p' 2d with c"(*> y) = p d+a c(p%, pv) ■ 

Remark: Note that we abuse our own notation here. The above definition of the family 
(C p )p>o does not correspond correctly to our sequence of conductivity functions (C n ) nG K 
defined in the introduction. To be precise: Given an arbitrary sequence (C n ) n( z^ in the 
sense of the introduction there might be no conductivity function C : Z d x Z d — ■> [0, oo) 
such that C p = C n for Vp = n. Nevertheless, we use C p in the sense above and C n in 
the sense of assumptions (Al) through (A5). This remark carries over to the definition on 
the family {Y n ) n& ^. We use the symbol Y n for the time continuous process corresponding 
to the conductivity function C n . That is, the family (Y p ) p>0 is determined by a single 
conductivity function C : Z d x Z d — > [0, oo) whereas the family (F n )„ gIN depends on the 
whole sequence (C™) ng M- Y n corresponds to the Dirichlet form (£ n , L 2 (n~ 1 'E d , n~ d p)) 
defined by 

(2-9) £ n (f,f) = \ E (f(x)-f(y)) 2 C n (x,y)n- 2d . 

x,yGn- 1 Z d 

Note that (Al) through (A4) are stable in the following sense. If one fixed conductivity 
function C satisfies (Al) ((A2), (A3) resp.) then the assumption holds true for the family 
(C p ) p with the constants independent of p. On the other hand, if (A4) is true for C the 
conclusion of (A4) holds for any C p with the same 7 and k 3 whenever r > 1 p~ 1 . 

Scaling as above implies a relation between the heat kernel pyp of Y p with respect to p~ d p 
and the heat kernel py of Y. Note that, by regarding the heat kernel of the scaled process 
with respect to p~ d /i, pvp{t,x,y) is not anymore the probability that the process starting 
in x is at time t in y but p~ d times this probability. One has 

(2.10) p Y p(t,x,y) = p d p Y (p a t,px,py). 

Let Y p,x be the process Y p with all jumps bigger than A removed. Y p,x corresponds to the 
Dirichlet form (£ P ' X , L 2 (p~ 1 Z d , p~ d p)) defined by 

^\fJ) = \ E (m-f(y)) 2 c p (x,y) P - 2d . 

\x-y\<\ 
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Finally, let V' 1 be the process associated to 

£ p /UJ) = \ E (/(*)- /GO) 



2 I i— d— a -2d 

ar-y p . 



x,y£p~ 1 'Z d 
\x-y\<l 



Let us finish this section with an overview over all processes which we have introduced so 
far: 



Y = (Y t 
Y p = (y/ 
Y n = (Y t n 
— ( 



Y p,x = (Yf' x 



V=(V t 
V p ' x = (V t p ' x 



: time-discrete Markov chain on Z d corresponding to the conductiv- 
ity function C. 

: time-continuous Markov chain on Z d with the same jumps as (X n ); 
its Dirichlet form is S. 

: scaled version of the process (Y t ); it corresponds to C p ; its state 
space is p _1 Z d ; its Dirichlet form is 8 P . 

: time-continuous process corresponding to C n ; its state space is 
rT Xr L d ] its Dirichlet form is S n defined in f |2.9j) . 

: limit of Y n for n —>■ oo; its state space is R, d ; corresponds to Dirich- 
let form (£,D(£)) defined in ( 01) . 

: equals the process (Y t p ) but with jumps greater than A removed; 
its Dirichlet form is S p,x . 

: rotational-invariant a-stable process on Z d ; its Dirichlet form is £ a . 

: scaled version of the process V with jumps bigger than A removed; 
its state space is p _1 Z d ; its Dirichlet form is S p ' x . 



3 Upper bounds for exit times and the heat kernel 

The aim of this section is to establish upper bounds on the heat kernel of the processes 
Y p,x independent of p > 1. This is done in section I3TTI The results are then applied in 
section 13.21 in order to establish Theorem 13.61 which is the key ingredient needed to show 
tightness of the family (Y n ) n( z]^. The techniques used in this section are borrowed from 
[CKS87] . [BL02bj . and \CKM\ . 
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3.1 Upper bounds on the heat kernel 



The main tools in this section are techniques worked out in [CK S87J. Lemma f3.ll is new 
and, together with Lemma I2.4[ we consider it important for the further development of 
anisotropic jump processes and Markov chains. 

Lemma 3.1 Assume (Al), (A2) and (A3). Then there exist c > depending only on N , 
k 2 and 02 > 1 depending on Ki, k 2 , No, d, a such that for all f G L 2 (Z a! , p), p > 0, A > 

S^U, f) < cE^{f, f) and m particular f) < c£ p (f, f) . 

Proof: Let {zf ,v ^) be the chains associated to C by (A3) now scaled on p _1 Z d . Note that 
(A2) and (A3) together imply \zj^ — zj x ' v ^\ < Q 2 \x — y\ with 6 2 = 6 2 (/ti, k 2 , N , d, a) for 
any chain in the sense of (A3), any pair (x,y) and any I. For notational convenience we 
assume the length of all chains to be equal to No. Then 

^'\fJ) = l E (m-f(y)) 2 \x-y\~ d ~ a p- 2d 

x,y£p- lT Z d 
\x-y\<X 

<n JL(f^-f)-M^)f\^-y\- d - a p- 2d 

x,y£p- 1 Z d 2=1 
\x-y\<\ 

No 

< n (k 2 )-i £ J2 (/(*&?) - f(4 x ' y) )) 2 c p (4 x J\ 4 x ' v) )p- 2d 

x.j/ep- 1 ^ 2=1 

\x-y\<\ 

< N*^)- 1 E p- 2d ^ { (f(ztf) - f{z\^)) 2 C{z\^\ z\^)} 

a;,2/Gp- 1 Z d 
\x-y\<\ 

<(Nof(K 2 y l 8^UJ). 

For the last inequality we use the fact that every term of the sum on the left appears at 
least once in the sum on the right hand side. By the second part of (A3) this happens at 
most Nq times. ■ 

Recall the following on-diagonal estimate for the truncated a-stable process on p~ lr L d given 
in [BL02bj : 
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Proposition 3.2 There exists c > independent of p > 1 such that 

Id-V" for f£ (0,1], 
" }„rt>l, 

i. e. in particular py P ,\ (t, x, y) < dr d l° c e t for all t > 0. 

These estimates lead almost directly to upper bounds for the heat kernel of Y p,x . 

Lemma 3.3 Assume (Al), (A2) and (A3). Let ©2 be the constant arising in the statement 
of Lemma \3.1\ Then there exists c > such that for all p > 1 and X > 2 



p YP ,x(t,x,y) < 



' ct -d/a f or t<l, 
ct -d/2 f or t>\. 



Proof: Applying Theorem 2.1 and Theorem 2.9 of [CKS87J to the assertions of Proposition 
13.21 one obtains 

2 r? a,d <c,(\\f\\%+£ P a \fJ)) 11/115", 
2 L t 4/d <c 2 (\\f\\ 2 L2 +£Z\fJ)) \\f\ff. 



Estimating S^' l (f, f) from above with the help of Lemma [3. II and applying afterwards the 
converse parts of Theorem 2.1 and Theorem 2.9 of [CKS87] the desired result follows. ■ 

In a similar fashion, one can use the scaling property ( 12.101) . Lemma 13.11 and the upper 
bounds on the rotational alpha-stable process on Z d (see Proposition 4.2 in [BL02b] ) to 
obtain directly upper bounds on the heat kernel of (Yf): 

Lemma 3.4 Assume (A2) and (A3). Then there exists a constant c > independent of 
p > such that 

p Y p{t,x,y) < ct~ d/a Vt > 0,^ x,y E p~ Xr L d . 
In fact, c only depends on d, a and the constants N and k 2 appearing in (A3). 

Off-diagonal estimates on the heat kernels of the truncated process can be obtained by 
Davies' method, see for example §3 in [CKS87J. 

Lemma 3.5 Assume (Al), (A2) and (A3). For any A > 02 there exists c > such that 
for all p>l,x,ye p~ Xr L d and t G (0, 1] 

p YP: x(t,x,y) < ct- d/a e-\ x - y \. 
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Proof: Fix x,y G p x 7^ y. Define 

r(/, /)(*)= E W*) - f(y)) 2 c p {x,y) P - 



yep- 1 ™ 4 



\x-y\<\ 

A 2 (^)=max{||e-^r(e^,e^)|| oo ,||e^r(e^,e^)|| oo }, 
£(t, x, y) = sup { \^{x) - i/>(y) | - tk^) 2 : A(^) < 00} . 

By Corollary 3.28 of [CKS87] and Lemma we get for t G (0, 1] 

p YP ,x{t,x,y) < ct - d /« e - E{2t > x > y) < ct-d/ae-WW-Tpm+ztAW 2 

for every ip with A(^) < 00. Fix x, y. Take ■?/>(£) = (j§5||,<0- Then t/>(x) — i/i(y) — \x — y\. 
Also we have \if>(£) - ?p(()\ < |f - CI- Now |e s - 1| < Ae A \s\ for \s\ < A. By (A2) we get 

e T2m T ( e ±^ je ±^ m = e T2m j2 (e**® - e^) V(c,0p _d 

l£-CI<A 

= ^ ( e ±wo-^))_i) 2 c"(c,0p~ d 
ie-ci<A 



< Kl AV A ]T |e-C| 2 "^ Q p^<c 4 A 4 -V A 



l?-CI<A 



where c 4 > is independent of x and y. 



3.2 Tightness 

We use the upper bounds established in the previous section to estimate exit times. For 
this we proceed as in [CK03] or [BL02b] and use the truncated process Y p ' together with 
a standard perturbation argument. 

Theorem 3.6 For any a > 0, b G (0, 1) there exists a constant 7 = 7(0, b, Ki, K2, N Q , d, a) 
such that for all R> 1 

P x (r(B 1 (x, aR); Y) < -fR a ) < b Wx G Z d and, equivalently, 
F x (r(B 1 (x,aR);Y p ) < jR a ) <b Vi G p~ Xr L d ,Mp > 1. 
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Proof: Fix A > G2. First, note that there exists by Lemma [3.51 a constant c\ > such 
that for all t G [1/2, 1] 

p YP ,x(t,x,y) < cie~ lx ~ yl . 
Therefore we estimate for t G [1/2, 1] and r > 

(3.1) F X (\Y P ' X -x\>r) = PYeAt,x,y)p~ d <c 1 ^ e^'V < c 2 e~ r ' 2 . 

\ x ~ y\~> r \ x ~ y\> r 

With the help of (13.11) we can now estimate the probability that Y p,x has left B p (x,r) 
before t < 1/2. Introduce the optional time T r := inf [t : \Y t p,x — Y P,X \ > r}, i.e. the first 
time the process has left the r-ball around its starting point. Then for all t < 1/2 we get 
by (13 .ip and the strong Markov property: 

P x (sup \Y P ' X - x\ > r) = F x (T r < t; \Y P ' X - x\ > J) + F x (T r < t; |Yf A - Y p ' x \ < J) 
s <t 2 2 

<F X (\Y P ' X -x\>-) + / F x (\Y 1 p ' x -Yf' x \> 7 -;T r eds) 
2 Jo 2 

< c 2 e- r / 4 + j t E x {F Y ^(\Y p l x s - Y p ' x \ > r ~;T r G ds)) 

< c 2 e~ r/4 + c 2 e- r/4 F x (T r <t)< c z e~ r/ \ 

Hence we have shown for all t G [0, 1/2] that 

(3.2) P x (sup \Y S P ' X -x\ > r) < c 3 e" r/4 . 

s<t 

We will now pass over to the process (Y t p ) by handling its large jumps as perturbations of 
(Yf ,x ) with standard techniques of perturbation theory for semigroups, see |Kat76j for the 
general case and |Lev72] for Markov semigroups. Let s^^ p ^ resp. J5fw be the generator 
of (Y t p,x ) resp. (Yf) and Q\ P,X ^ resp. P/ be the corresponding semigroups. Then 

jgfto = + gs^ x) 

where 

x,y£p- lr E d 
\x-y\>\ 
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s. 



This is a bounded operator in L°°(p 1 Z d ) since for bounded / by (A2) 

\-^ x) f(x)\< J2 \f\y)-f^)\c p (x J y)p- d <2K 1 \\f\\ 00 £ k-yr d_Q p- d 

x,y£p- 1 Z d x,y<=p- 1 Z d 
\x-y\>\ \x-y\>\ 

<C4||/IL. 

Define Sb(t) = <?S P ' A) and for k > inductively 

W*) = f S k (t-s)^S (s)d 
Jo 

Then it is immediate by the contraction property of Q^'^ that Sfc(£) is again bounded in 
L°° with operator norm 

ll&wiL. s 

Clearly X)fcLo^fc(0 ^ s well-defined and equals the perturbed semigroup Pj: p \ We have 
||pW -Q^' A) ^ = || £~ Q Skit)^ x < c A te CAt and can therefore estimate for t E [1/2, 1] 

F x (\Yf -x\>r)= P t (p) t BP(X:rr < Q { t pA) t BP(X:rr + c 4 te Cit = F X (\Y P ' X - x\ > r) + c 4 te Cit 
< c 5 e~ r/2 + c 5 t. 

We now proceed as above. Set T r := inf [t : \Y P — Y P \ > r}. Then for all t < 1 we get 
again by the strong Markov property: 

P*(sup \Yf -x\>r)= F x (T r < t; \Y P - x\ > + F x (T r < t; \Y P - Y p \ < 



s<t 



< F X (\Y P - x\ > T -) + j\ x {\Y p - Yf\ >^,T r e ds) 

< C5 e-r/± + C5t + J w(v Y '\\Y{_ a -Y*\> r ^T r eds)) 

< c 5 e- r/4 + c 5 t + (c 5 e" r/4 + c 5 t)F x (T r <t)< c 6 e- r/4 + c 6 t. 



This translates by scaling into the following estimate for the process Y: 

F x ( sup \Y S - x\ > pr) < c 6 e~ r/4 + c 6 t. 

s<p a t 

Given a and b we choose r > and t < 1/2 such that the left hand side is smaller than b 
and in addition a/r > 1. Now setting p = aR/r proves our claim with 7 = a a t/r a . ■ 



For later purposes, choose 7 = 7(1, 1/2), i.e. 

(3.3) F x (r(B 1 (x,r);Y)<^r a )<^. 
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4 Hitting time estimates and the regularity of the 
heat kernel 



In this section we derive an equicontinuity result for the heat kernels of the processes Y p . 
In our application it is essential that the constants appearing do not depend on the scaling 
parameter p > 1. Again, our presentation uses results from |BL02bj and |CK03j . Another 
option would be to adopt methods of |Kom95j . 

First observe the following Levy system identity, cf. [CK03J: 

Lemma 4.1 Let f: R + xp~ 1 Z d xp _1 Z d — ► R + be a bounded measurable function vanishing 
on the diagonal, i.e. f(t,x,x) = for all x G p~ lr L d . Then for all x G p~ lr L d and 
predictable stopping times T we have 

E" Yf_,Yf) 

Let W p = (Wf)t be the space-time process on R + x p _1 Z d associated to Y p , i.e. W t = 
(Ut,Y t p ) where Ut = Uq + t is a deterministic process. We call a measurable function 
u: R + x p _1 Z d — > R space-time harmonic or caloric on an open set Q C R + x p _1 Z d if for 
all open relative compact sets Q' C Q, (t, x) G fl' 

u(t,x) = E^(u(W^ wn )). 
Important examples for space-time harmonic functions are given by the heat kernel of Y p : 

Lemma 4.2 Let to > 0, y G p _1 Z d . Then the function u(t,x) = pyp(to ~ t,x,y) is 
space-time harmonic in [0,to) x p~ lr L d . 

The proof is exactly the same as in Lemma 4.5 of [CK03] . 

Define Q p (t,x,r) := [t,t + ^r a ] x B p (x,r) where 7 is chosen such that (13.31) holds. We 
have the following estimate on the probability of hitting relatively large sets before exiting 
Q p (0,x,r): 

Lemma 4.3 Assume (Al), (A2), (A3) and (A4). Then there exists a constant c(K^,d, a) > 
such that for any r > ©ip~\ x G p _1 Z d and compact set A C Q p (0,x, r) with A® \i p '(A) > 
l\®H p (Q p (0,x,r)) 

F^(a{A;W t p ) < T(Q p (0,x,r);W t p )) > c. 



E 1 



f(s,Y s p ,y)C p (Y s p ,y)p- d )ds 
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Remark: More general, with the proof below and Gi, 7 as in Lemma [231 we can get lower 
bounds on the probability of hitting sets A C Q p (0,x, r) with A ® p p {A) > 2(1 — 7)A ® 
p(Q p (0, 2, r)) before exiting Q p (0,x,r) for r > ©ip -1 only with assumptions (Al)-(A3). 
Unluckily, for technical reasons we need (A4) to prove our equicontinuity result. 

Proof: Set r r = T(Q p (0,x,r);W t "), o A = a(A;W t ") and T = a A A r r . Without loss of 
generality we may assume F i0 ' x \a A < t t ) < \. For each s E [0, 00) let A s denote the 
projection of A on {s} x p~ lr E d and let for y G p~ lr E d 

N(y) = {ze p-'Z d : C(y, z) < k 3 \y - z^} . 



Choosing /(s,£,C) = t QP{0tXtr)xAs \ {{yty) . yep -i zd} (^ () in the Levy system formula implies 



E(o.») 



> E (0 ' s) 



E CP{Y ^ z )p~ d ) ds 



z&A s \N(Yf) 



> k 3 E^ x) 



E w 



\—d—a ~d \ 7 

z\ p as 



> 2- d - a K^ x) 



z£A s \N(Y s p ) 
■T 







E r~ d - a p- d \ ds 



z£A s \N(Y s p ) 

> 2- d - a K 3 r- d - a E {0 > x) \ [ p p {A s \ N{Y p )) ds 



> 2- d ~ a K 3 r~ d - a 'E^ x) 
1 



p p (A s \ N(Y S P )) ds;cr A AT r > ^r a 



> — ■ 2 

36 



-d— a . 



K 3 r- d ~ a (X^ p)(Q p (0,x,r))F ( -°' x) (a A Ar r > ^r a ) 
>c(K 3 ,d,a)F^(a A AT r >l^r a ). 



Here we have used (A® p)(Q p (0, x, r)) x r +a . For the second last step note that for every 



path with a A A r r > f tt 



-\®p"(Q p (0, x, r)) < A <g> pL" {A) 



1 



< j p p (A s \N(Y s p ))ds + J p p (A s nN(Y p ))ds + -A® p, p (Q p (0,x,r)). 
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Now we get by (A4) and because of r > 3 p d 
J » P (A S H N(Yf))ds < J p(B p {Y p ,r) n iV(n p ))ds <^J fi p (B p (Yf,r))ds 



Hence 



^y(^(x,r)) = ^\®p p (Q p (0,x,r)). 



-"IT 1 

\ N(Y p ))ds > —X (g> p p (Q p (0, x,r)). 



By our choice of 7 we obtain 

F i0 ' x \r r < §7r a ) < F x (r{B p (x,r)]Y p ) < ^r a ) < -. 

Finally we estimate 

F {0 ' x \a A A r r > §77-") = 1 - F^°' x \a A A r r < \^r a ) 

> 1 _ F^ x \a A < T r ) - P (M (tv < ^r a ) > ~. 



We also need the following upper bound on the probability of exiting a ball by large jumps. 



Lemma 4.4 Assume (Al), (A2) and (A3). Let 0i > 1 be the constant of Lemma \2.4 



Then there exists a constant c(k\, K2, N Q , d, a) > such that for all p > 1, s > 2r > Q\p 1 , 
(t,x) E [0,oo) x p- lr L d 



p(t,«) (\V P 



(QP(t,x,r);WP) 



£ Q(t,x,s)) < c— . 



Proof: Let r = r(B p (x, r); y). Observe that |£ - C| > §|z-C| for £ e B p (x,r) and 
£ G B p (x, s) c . Since the space-time process moves continuously in time, it can only exit 
Q p (t,x,s) and Q p (t,x,r) simultaneously by jumping in space. Using this fact together 
with the Levy system identity for (Y t p ) and (A2) one obtains 

V M ( W «QP<&,r);WP) i Q P (t,X>*)) = ^ ( Y r t B P {x i s) ] T< J^) < F X (Y p £ B p (x,s)) 

J2 C p {Y p ,y)p- d )dt 

\y-x\>s 

d-a p - d 



E X [ ^2 ^-{Y t p _£BP(x,r),Y t p £BP{x,s)} 

t<T 



= E x 









<wyj K i \ Y t - y\~ d ~ a p~ d dt\ < ^i2 d+a E^r \ x ~y\ 

\v-x\>s \y-x\>s 
< ds _a E*(r). 
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Therefore it remains to estimate E s (t). Let K3 and 7 be the constants of Lemma 12.41 
For each f G p- lr L d define A{g) = {( G p^Z d : C"(f,C) > «s If - Cl~ d " a }• We get by 
applying the Levy system identity in the above fashion for r > ©ip -1 



1 = F*(Yf i B p (x,r)) = w[YhYf^B {x ,r),Yfmx. 

t<T 



'■)} 



JO 



E' 1 



> E* 



> k 3 E x 



I j/ — a;|>7- 



^ C(Yf,y)p- d )dt 



\y-Y t p \>2r 
y&A(Y t ") 



2r<|?/-y t p |<6r7- 1 / d 

veA(Y t p ) 



> fc 3 E a 



£ CO(Yf,y)p- d )dt 

\y-Y t "\>2r 

b-F/|>2r 
y&A(Yf) 



> K 3 7- d ~ a -? 1+d/a E x / p"V(^(^ P ) H B p (Yf, 6n~ 1/d ) \ B P (Y P , 2r))dt 

-Jo 

> d(/C3, tZ, a)E* [ / p- d p(B p (Y t p , 2r))dt\ > c 2 (k 3 , d, a)r- a W(r). 



Here the second last inequality is due to Lemma 12.41 since there are at least 2p(B p (Yf, 2r)) 
elements of A(Y P ) in B P (Y P , 6r~/- l / d ). ■ 



Proposition 4.5 Assume (Al)-(AJ^). Then there exist constants c > and (3 G (0,1) 
depending only on the constants appearing in (A1)-(A4) such that for all R > ©ip -1 , 
q: [0,716-R] x p~ lr L d — > R bounded and space-time harmonic in Q p (0,Xq, 16R) the following 
a-priori continuity estimate holds: 

\q(s,x) -q(t,y)\ < c WqW^ FT? ( \s - t\ 1/a + \x-y\) 

for all (s, x), (t, y) G Q p {0, x , R) with \x — y\ > ©ip -1 . 

Moreover, if \x — y\ < ©lp" 1 we have 

\q(s, x) - q(t, y)\ < c || g |L JT'( \s - t\ l/a + B 1 p- 1 f . 

Proof: The proof can be found in [BK06] or |CK03j for our case and in [BK05] or [HK06J 
for the "elliptic" case. 

We may assume HgH^ = 1/2. Fix (t,x) G Q p (0,x ,R). Let £ G (0,1/2), 77 > 0, Q k = 
Q p (t,x ,£ k R), r k = r(Q k ; Z p ) and define for k G IN 

m k = inf q(z), M k = sup q{z). 
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We show that it is possible to choose constants £, ( independent of R > 0, x, x G 
Q p (0, x , R) and q such that 

(4.1) M k -m k <( k 

for all k > with > 0ip _1 where Gi the constant in (A4). The restriction that x and y 
cannot be arbitrarily close is natural since hitting time estimates in the sense of Lemma |4~31 
may not hold for small r. Just consider the conductivities C(x, y) = \x — y\ a l{| x _ y |>^ }. 
Then the process {W l ) has only jumps with length bigger than Rq and therefore can't hit 
sets A C B(xo, r) for r < Ro starting in xq unless xq G A. 

Trivially, (14. Xp holds for k = 0. Now assume that this equation holds already for alH < k 
while still ^ k+1 R > Q^ 1 . We set 

A k = [z G Q k+ i : q(z) < }. 

Without loss of generality we might assume \<g> n(A k ) / \<g> /i(Q k+ i) > 1/2. Else we just look 
at 1/2 — q instead of q. We choose a compact set A' k C A k with A <g> n(A' k )/\® fi(Q k+ i) > 
1/3 and define T k = T(A' k ;W p ). Observe that u is harmonic in Q k C Q p (t,x,R) C 
Q p (0, Xq, 16R). Therefore we get for arbitrary Zi,z 2 G Q k +i 

q(zt) - q(z 2 ) = E 21 [q(W P kATk+i )] - q(z 2 ) 

= E 21 [q{W P k )-q{z 2 )-T k <r k+1 ] 
+ E 21 [q(W Tk+1 ) - q(z 2 ); T k > r k+1 and W p k+i G Q k ] 

k 

+ £ E 21 [q(W Tk+1 ) - 9(22); T fe > r k+1 and G Q fe _ 4 \ Q k - l+l ] 

i=l 

+ E 21 [g(W Tft+1 ) - g(^); T fc > r k+1 and £ Q ] . 

Here, the first term can be estimated by \{M k — m k )V zl (T k > r k+ i) while the second term 
is bounded from above by (M k - m k )F Zl (T k > r k+1 ) = (M k - m k )(l - P 2l (T fc < r k+1 )). 
Moreover Lemma [4.31 implies the existence of c% > such that P 2l (Tfc < r k+ \) > c±. For 
the remaining terms note that by Lemma 14.41 there exists c 2 > with 

V Zl (W T p k+i iQi)<c 2 ^ k+1 -^ a 
Summing up the terms above we this get for £ a < (/2 

, k+l 
q{zi) - q(z 2 ) < {M k - m k )(l - -P 21 (T fc < r fc+1 )) + ^(M fc _, - m fc _,)P 21 (^ +1 £ Q k - i+ i) 
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<(i-I Cl )c fe -^ic fc +c 2 ^c^r 

i=i 

< (1 - \ciX k - \dC, k + c 2 ( k f; = (1 - ^ Cl )C* - \ Cl ( k + 2c 2 £( k . 

1=1 ' 

Estimate (14. ip follows from the equation above by taking ( = 1 — ci/4 and £ = ~ A (I) 1 a A 

Cl 

8c 2 ' 

To derive Holder continuity let Zi = (s^Xi) G Q p (0, x , R), i — 1,2 with ^ 7^ z 2 and 
Si < s 2 . Assume |xi — a^l > @iP _1 and take k maximal such that z 2 G B p (zi, Rl; k ). Then 

\si - s 2 \ 1/a + \ Xl - x 2 \ < (f' a + l)Re, eR > ©iP~\ \q{zi) - q(z 2 )\ < ( k . 

Thus by optimality k is the smallest integer such that 

k > (logO'^logdsi - s 2 \ l/a + \ Xl - x 2 \) - log(7 1/a + l)R) - 1, 

and we get 

\q{zi) - q{z 2 )\ < C\t /a + l)- lo ^/ lo ^R- lo ^ lo ^{\ Sl - s 2 \ 1/a + \ Xl - x 2 |) logC/log «, 

i.e. the proposition holds with (3 = log£/log£ g (0,1) and c = C~ 1 (7 1 ^° + l)"' 3 - If 
on the other hand \x\ — x 2 \ < @i(7 1 / a + l)p _1 we take k maximal with R£ k > 0ip -1 
and z 2 G B p (z\ } Rt; k ). Then in particular \s± — s 2 | < 7 1 / a i?^ fc , and we get for k the 
inequalities 

k< (logO-^iogCei^-iogiz), 

k < (logO-^logd^ - s 2 \ 1/a ) - log(7 1/a i?)) 
Combining this we get with (3 as above 

\q{zi) - q{z 2 )\ < C\r P,a + l)R~ P ( \si - s 2 f /a + (Oxp- 1 )^) 



In particular, this implies regularity of the heat kernels 
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Theorem 4.6 There exist constants c > and (3 G (0,1) such that for all to G (0, oo), 
x ii x 2, V £ p~ Xr L d and s 1; s% > [t , oo) For y G p _1 Z d we /ictue 

|pyp(si,aci,2/) -PKf(s2,a?2,3/)| 

< c t (d+/3)/Q ( | ai - s 2 \ 1/a + \ Xl - x 2 \ v (e^ 1 ))^. 

More general, for arbitrary yi,y 2 G p _1 Z d 

\PYp(suXi,yi) - PYp(s 2 ,x 2 ,y 2 )\ 

< ct- {d+p)/a {\ Sl - s 2 \ 1/a + \xi - x 2 \ V (Q 1P - 1 ) + \ Vl - y 2 \ V (Q 1P - 1 )Y. 

Proof: Fix to > 0. For an arbitrary T > to the function q(t, x) = p(T — t, x, y) is space- 
time harmonic on [0, To/2) x p~ lr L d by Lemma l4~2l as well as bounded by c\t Q d ^ a by Lemma 
031 Now take R = {T /32*/) 1 / a and s u s 2 G [0,7-R a ). Assume first Qip' 1 < |xi - x 2 \. If 
l^i — x 2 1 < R we get by Proposition 14.51 

\p Y p(T - s x ,xi,y) -p Y p(T - s 2 ,^ 2) y)| < c 2 (T /32j)-P /a c 1 t~ d/a (\si - s 2 \ 1/a + \xt - x 2 \ ) 

<c 3 t - (d+/3)/a (k 1 - S2 r /Q +ix 1 -x 2 |) /3 . 

In the other case \x± — x 2 \ > R we have (|si — s^ 1 ^" + |#i — x 2 \Y > c^J a and hence 

\p Y p(T -s u xi,y) -p Y p(T - s 2 ,^ 2 ,y)| < 2cit a!/a < 2cit ( d +^)/ Q tJ/ a 

< c 5 to (a!+/3)/Q ( |si - s 2 | 1/a + \xi - x 2 \ ) . 

In the same fashion we can deal with the case 0ip _1 > \xi — x 2 \. The other a-priori 
estimate asserted in the theorem now follows by symmetry of the heat kernels. ■ 



5 The central limit theorem 

The aim of this section is to provide a proof of Theorem 12.11 and Theorem 12.21 Although 
the limit process SC is a jump process the idea of the proof is very similar to the one in 
[SZ97J and [BK06J. For n G IN let Y n = (Y t n ) t be the time- continuous processes defined 
in the introduction and explained in section 12.41 Denote the corresponding semigroup by 
(P t )t and its kernel by p( n \t,x,y). Recall that the Dirichlet form corresponding to Y n 
is given by 

£ {n \f,f) = £ (/M - f(x)) 2 C n (x,y)n- 2d . 

x,y£n- 1 Z d 
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We denote the restriction of functions on R d to n _1 Z by R n . We also need to extend 
functions on the grid to continuous functions on R d in a (for our purpose) reasonable way: 
For x G rT lr L d set n (x) := n[xj,Xj + V n ) an d &n = U {0. n (x)}. The sets of n 

form a partition of R d . Recall that, for a point x G R d we denote by [x] n the element of 
n~ Xr L d satisfying ([x] n )j = n _1 \nxi\ for all i = 1, . . . , d. For /: n~ lr L d —>■ R we denote by 
_E n / a Lipschitz-continuous function .E n / : R d — > R satisfying: 

a) (x) = /(z) Vzerr'Z", 

6) /(^) < E n f(x) < f(x) Vx G R d , where /(a;) = min /; /(x) = max / , 

c) \VE n f(x)\ < cnmax{|/(£) - f(rj)\ : £, v G O n (i)flC 1 Z li } Vn G !N,Vx G R d . 

The precise choice of the function i? n / is not important for our approach as long as E n is 
a linear operator. 

Let us emphasize that in the following result we adopt the notion P x for the probability 
of a Markov process starting in x. Any starting point, together with a stochastic process 
with cadlag paths, corresponds to a probability measure on D([0, oo), R d ). For a Markov 
process X = (X t ) starting in x we refer to this probability measure as "the law of X under 

The following theorem is a reformulation of Theorem 12.11 

Theorem 5.1 Let (C n ) n be a sequence of conductivity functions satisfying (Al) through 
(A4) and (x n ) n a sequence of points x n G n~ lr E d with x n — > x G R d forn —>■ oo. T/ien eac/i 
subsequence (n') of (n) has a subsequence (n") with the following properties: 

1. For any f G C c (R d ) the continuous functions (^E n rrP^ n ^R n »f) converge uniformly 
on compact sets for n" — > oo. The limit defines a family of linear operators (Pt)t>o 
which extends to the semigroup on C(R d ) of a symmetric strong Markov process 3£ . 

2. For any t > the laws of (Y t n ")te{o,t ] under P x ™" converge weakly to the law of 
(&t)te[o,ta]- 

Once these assertions are proved it remains to show that does not depend on the 
choice of (n'). The proof of Theorem 15. II makes use of the following sufficient condition for 
tightness: 

Theorem 5.2 (|A ld78| ) Let (Y n ) be a sequence of stochastic processes with cadlag paths, 
t > 0. Assume: 
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1. For all sequences a n of random variables with values in [0,to] such that a n is a 
stopping time with respect to the natural filtration of Y n , sequences 5 n > with 
lim 5 n = and r\ > it holds that 

(5-1) limF(\Y: n+Sn -YZ\>ri)=0. 

2. Either (Yq 1 ) and max^o^] \ Y t n — Y£_\ are tight or Y t n is tight for every t G [0, to]- 



Then the laws of (Y n ) n are tight in D([Q,t ], R d ). 



Proof of Theorem 15.11 and Theorem I2.lt Let (n') be a subsequence of (n). Fix 
countable dense subsets (sj) of [0, oo) and (fj) of C c (R d ). 

By := E n PtRn we define a positivity-preserving contraction semigroup {Q^)t on 
the Banach space C(R ). Now, Theorem 14.61 yields that for all i, j the family of functions 
(Qs? fj)neti is equicontinuous. In fact, we have for x, y G n~ lr L d with |x — y| > n -1 ©! 



< cs 7(^)/a A(supp(/ . )) || / .|| oo | a; _ y |^ 

where c > 0, Oi > and /3 G (0, 1) are independent of n and stem from Proposition 14.51 
The construction of the extension operators E n implies for all x, y G R d with \x—y\ > n. _1 03 

(5.2) IgW/iOO - < 2« ( -^ a A(supp(/ i )) ||/, IL |x - 

In a similar fashion one establishes for \x — y\ < n~ 1 Qi 

(5-3) \Q<£ ] m - Q { : } fM\ < 2ce lS r(^A(supp(/,)) \\f<\\__ n-f* 



J lloo 



Furthermore (Qs? fj) is equibounded. The Theorem of Arzela-Ascoli and the passage to 
a diagonal sequence give us therefore a subsequence (n") of (n') such that for all i, j 
the sequence Qs? /j converges uniformly on compact sets for n" — > oo. Denote this limit 
function by P Si fj- We use an £ /3-argument to extend it for all positive times: Let t G (0, oo) 
and take a subsequence (i r ) of (z) such that — > t for il — > oo and Sj/ > t/2. Then 

\Ql n " ] m - Qi m,,) m\ < \QPm - Qfpm\ + \Q%' ] m - Q^" ] m\ 

+ \Qt" ) f^)-Qt" ) fM- 
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The second term on the right hand side converges uniformly on compact sets to for 
n",m" — > oo. The other two terms can be handled again by Theorem 14.61 

\P t in] R n (fj)(x)-P£ } Rn(mx)\ < \p in) (t,x,z)-p^(s l ,x,z)\\f J (z)\n- d 

< W^/^supp^)) WfAL (I* - t\ 1/a + e^f. 

The right hand side clearly converges to for n — > oo. Hence the limit P t fj exists uniformly 

II in") II 

on compact sets for all t G [0, oo). Finally by \\Ql fj\\ < \\fj\\ and because (fj) is dense 
in C(K d ) in the topology of uniform convergence on compact sets we have established the 
desired convergence result for all / G C(¥L d ). 

It follows from the corresponding properties of the Q^ 1 that Pt is a positivity-preserving 
contraction semigroup on C(R d ) which is hence associated to a symmetric strong Markov 
process S£ on R d . 

Fix to > and x G R d . We want to apply Theorem 15.21 Take an arbitrary sequence of 
stopping times r n G [0, to], a sequence (S n ) of reals converging to and a > 0. By Theorem 
13.61 for each choice b G (0, 1) there exist a constant 7(0, b) with 

(5.4) P x " (r(B(x n , a); Y {n) ) < 7(0, b)) < b 

for all n G IN. Therefore, for all n large enough such that 5 n < 7(0, b), 

^(\y}:k - Y£>\ >a)= P-(|yW - Y ^\ > a) 

<P x "{r(B(x n ,a);Y^)<5 n ) 

< (r(B(x n , a); F (n) ) < 7 (a, 6)) < 6 

where the first equality follows by the strong Markov property. (15.11) follows immediately. 
Moreover, x n — > x implies the tightness of the starting distributions while (15. 4p implies the 
tightness of max te [ 0A) ] \ Y} n ^ — Y t _ |, both under F Xn . The tightness of the laws of {Y^) 
under V Xn follows. 

Finally we prove the asserted weak convergence for n" — » 00 by showing that the finite 
dimensional distributions of the limit probability Q on D([0, t ], R d ) of a weakly convergent 
subsequence (n'") are independent of the actual subsequence. For g G C c (R d ), t G [0,t ] 

[ g(uj t )d£>(uj) = lim W^"{R nlll g){X { t n '" ) ) = lim Py"\R n ,„g){x n ,„) 

J n'"—*oo n"'-^oo 

= Pt.g(x) , 

in'") 

where the last equality follows from the equicontinuity of the family P t g. Therefore the 
one-dimensional distributions are independent of (n'"). More generally let < si < . . . < 
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Sk < an d gx, ■ ■ ■ , gu G C c (H d ). Then by the time-homogeneity of our Markov chains 
/ gi {uj sl ) • . . . • g k (u,. k ) dQ{u>) = lim W n "' {g x {X^P) • . . . • g k {X Sk )^)) 

= (Pax (9lP S2 - Sl (• • • Psk-s^gk))) {x) ■ 

Here the last equality is again due to the equicontinuity. Hence the k-dimensional distri- 
butions of Q are independent of the choice of the subsequence (nf") and are determined 
by the semigroup (-P*). Therefore we have weak convergence along (n"). In particular, the 
stochastic process corresponding to Q has the same finite-dimensional distributions as 2£ 
starting in x. ■ 

We proceed to the proof of Theorem 12.21 

Proof of Theorem 12.21 Let (n') be any subsequence of (n). Let be a strong Markov 
process - possibly depending on the choice of (n') - and {n") be a subsequence of {n') such 
that the assertions of Theorem 15.11 hold true. We aim to show that SC does not depend 
on the choice of (n'). It suffices to show that the limiting process S£ corresponds to the 
Dirichlet form (II. ip . This is the case if 

(5.5) £(U x f,g) = (f,g)-\(U x f,g) 

for any f,ge C c °°(R d ) where U x f{x) = J °° e - xt {P t f){x)dt, A > 0. Note that, at this stage, 
U\ does depend on the choice of {n'). Equality (15.51) implies 

£(U x f, g) + X(U x f, g) = £(G x f, g) + X(G x f, g) for any /, g G C?(R d ) , 

where G x and (S',D(S')) are independent of {n'). G x is then the L 2 -resolvent of and 
we are done. Note that Theorem 12.91 implies 

D{£) = H a/2 (R d ). 

We prove (15. 5p by approximating each term by its discrete analog. On the discrete level 

(5.6) £^{U^R n {f),R n {g)) = (R n (f),R n (g)) - A(C/^ R n (f) , R n {g)) 

where U^h(x) = J °° e~ xt (pj n 'fi) (x) dt, A > 0, denotes the resolvent of (.Z^ ). 

Therefore fix A > 0, /, g G C^°(R d ) and abbreviate f n = R n (f), g n = R n (g)- Then 
/„, g n G L 2 (n _1 Z d , n _1 /i) with ||/ n || + ||g n || < c for all n. Recalling the definition of section 
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Oone sees that ^ ien -iz" fn{x)t Qn {x) converges in L 2 (R d ) to / and \{f n ,9n) ~ (/,</) | 
converges to zero for n — > oo. 

Now by the compactness of the support of / and Theorem 14. 61 we get equicontinuity for the 
family (E n U^ R n f) n &¥i analogous to (I5.2p resp. (I5.3p . Together with \[x] n — x\ < \fdn~ x 
we get 

\E n U^f n (x)-UPf n ([x} n )\ <cn- p 

for all x G R d . In particular, the functions x i— > U^ Vn'(Mn') on R d converge along 
the subsequence {n") uniformly on compact sets to U\f . Taking into account that g is 
compactly supported we get by dominated convergence 

(5.7) | (Uf Vn", g n ») - (U x f, g)\^0 for n" - oo. 

Therefore the right-hand side of ( 15. 6j) converges against the right-hand side of ( 15. 5ft for the 
subsequence n" — >■ oo. 

It remains to determine the limit of the left-hand side of (15.61) for n" — > oo. We do this in 
several steps. 

Step 1: U x f E H a ' 2 {R d ). 

This result probably follows from standard arguments of approximation theory. For the 
sake of completeness we give a detailed proof. First note that U^' f n and E n U^ f n form 
bounded sequences in L 2 (R d ). Set F n = E n U^f n . Then we aim to prove ||-Fn||,HW 2 (R d ) < c 
with c > independent of n. Define V n = {z G R d : |z|oo < 2n~ 1 }. Moreover, let 
. . . , z& be the corners of n (O). We write 

i._,.^a — rf ^C = 2^ // 

+ £ // (^(e + ^-F„(0) 2 |r ? r d - a rfr/^=:(/i) + (/2) 

xen " 1Zd Q n (x)x(R'«\Vn) 

Let us first look at (Ji). For x G n _1 Z d , £ G 77 G V^, by Taylor's formula 

(F n (e + ^)-Fn(0)' a \- (F n (J)-F n (y)) 2 

< c n > 



d+a — ^ u '" |„|d+a-2 

ireQ n (x)nn" 



yeQ„(x)+v n nn" 1 z d 

Furthermore, 



Qn(i)xK | ?? |<2v / dn- 1 
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where c\ > only depends on a and d. Since \y — x\ < A\fdn 1 the first sum in (15.81) can 
be bounded from above by 

c * n ~ 2d E E {F n {x)-F n {y)) 2 \x-y\- d - a 

yeQ„(x)+F„nn- 1 z d 

x^y 

<c 3 n- 2d {Fn(x + h)-F n (x)) 2 \h\- d ~ a . 

x,h^n- 1 Z d 
0<|/i| oo <6n- 1 

As expected, (ii) tends to zero for large n. In order to tackle (J 2 ) note that, for all 
h, x e n _1 Z d , 7] e H n (x + h),£e 0. n (x) 

(FM - F n (£)) 2 < max (F n (x + h + z^) - F n (x + z[ n) )) 2 

i,j = l,...,2 d ' 

< J2 {Fn(x + h + zf ] ) - F n (x + zf)) 2 . 

hi 

For any Q. n (h) C R d \ V n , h G Z d , and any r? e 0. n (h) and i, j = 1, . . . , 2 d 

\V\ > Moo > \\h + - L > ^ + *< n) - ^| > c*|/i + zt ] - zf\ 

where C4 depends only on the dimension. Keeping in mind that the volume of H n (x) x 
0. n (x + h) is rT 2d we estimate (I 2 ) in (15 .8p from above by 

c2 +d n~ 2d jr (F n (x + h + zt ] ) - F n (x + zf)) 2 \h + zf - zf^ 

M=l x,h£n- 1 'Z d 

, . (n) i n ) _Lr\ 

<c 5 n- 2d {F n {x + h) - F n {x)) 2 \h\- d ~ a . 

x,h£n- 1 T, d 

Hence we obtain by (A3) and ( 15. 6ft 

R d xR d 

<c 5 n" 2d £ {U^ f n (x + h) - f n (x)) 2 \h\- d - a 

x,hen~ 1 Z d 

<c 6 S^(U^f n ,U^f n ) =c 6 (U^f n ,f n ) - c 6 X\\uP f n \\ 2 . 
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The right-hand side is bounded in n. We conclude that E n U x f n is a bounded sequence 
in the Sobolev space H a / 2 (R d ). 

In our situation this means that there is a subsequence (n'") of (n") such that E n irtU^ ^f n >» 
converges weakly in H a / 2 (ja d ) and strongly in L 2 (K) for K C R d compact to an element 
F G # a / 2 (R d ) for n'" -> oo. Since E n »U { x n " ] ' f n „ -> C/ A / pointwise, C/ A / = F almost 
everywhere. In particular C/ A / G iJ Q / 2 (R d ). 

Step 2: Setting for r G (0, 1) £ r (/, g) = \ J J (f(y) - f(x))(g(y) - g(x))k(x, y) dx dy one 

|x-y|>r 

observes 

(5.9) \S r (U x f,g)-S(U x f,g)\^0 for r - 0. 



Step 3: In analogy to Step 2 set 

-2d 

2 

|s-i/|>r 



x,y£n- 1 '* d 



Then for any r > 



(5.10) 



£("")([/f <?n") - £r(Uxf,g)\ - for n" -> oo 



Assertion ( 15. 10D is easily established by estimating the difference for large enough n" from 
above by 

(^f'V([y]n») - f/f'V (Mn'O) ( 5 (Mn») " <7 (Mn»)) C< , fa/]„«) 



\x-y\>r/2 



(U x f(y) - U x f(x))(g(y)-g(x))k(x,y) 



dx dy 



This term tends to zero for n" — > since g{[x] n ») — > #(:r) and V(Mn") "~ > U x f(x) 
uniformly on compacts for n" — > oo. Using (A5) we see that the integrand converges 
uniformly to 0. (15.101) follows by dominated convergence. 

Step 4: 

(5.11) sup \£^ (Uf g n >) - £f ] (U^fn^gw) |-0 for r - 0. 



Recall that both, / and g have compact support, i.e. the number of elements in supp(g) fl 
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n 1 Z d is of order n 2 d. Using Cauchy-Schwartz and (A2) we obtain 



\x-y\<r 
x,y£n- 1 Z d 



E (U^fnix) ' U^Uy))(g n (x)-g n (y))C n (x,y)n 

<( E (u^f n { x )-u^Uy)) 2 c n {x,y) 

\x~y\<r 
x,y£n~ 1 Z d 

x E {9n(x)-g n (y)) 2 C n (x,y)n- 2d 

\x-y\<r 



n 2d 



x,yGn~ 1 Z d 



< 



c\S^(U^f n ,uPf n )\\S r (g n ,g n )\. 



Let us look at the above estimate for n = n". Since |£( n ")([/j n ^/ n «, [/■[" Vn")| and 
\£{9n", g n ") | are bounded uniformly in n" 



sup |£("")(C/f Vn",^" } /n»)| |£r(<7n»,<7n")| - for r ^ . 



Step 4 is completed. 

Finally, combining (15. 9p . (I5.10p and (15. lip by a standard chaining argument proves (15. 5p . 



6 Approximation of jump processes by Markov chains 

Here, we prove Theorem 12.31 

Proof of Theorem 12.31 Obviously, by the symmetry of k (Al) holds. The upper bounds 
on k imply for x, y G n _1 Z d , \x — y\ 00 > 2n~ 1 : 

C n (x,y)<K 5 n 2d J \Z-(\- d -«dZd( 

|x-g| oc <n- 1 /2 
|»-Cloo<n _1 /2 

< K 5 {4dY d+a V 2 n 2d [ \x-y\- d - a d£d( < ^{Ad) {d+a ^ 2 \x - y\- d ~ a . 



|y-CL<n-V2 



since |f - C| > |f - CL > \ x ~ vloo ~ n 1 > \ \ x ~ vloo > \ d 1/2 \ x ~ v\- In the same wa Y 
one shows C n (x,y) > Ki(Ad)~( d+a ^ 2 \x — y\ ~ d ~ a for \x — y\ OQ > 2n _1 . Therefore, (A2), 
(A3) and (A4) are satisfied. 
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Let k n (x,y) = C n ([x] n , [y] n ). Fix a compact rectangular subset Q = Yli=i[ a iih]i a i < k of 
R d x R d \ diag and define Q e = n^=i[ a « ~ £ ,h+ e] for e > 0. In particular, there is c\ > 
with A(fi e \ fl) < Cie for e < 1. Then there exists n Q > such that dist(fl, diag) > rig 1 . 
For all n > Ano and e < (4no) _1 the functions k n are uniformly bounded from above on 
fl £ by C2 > 0. By the Theorem of Lusin, for each e > there exists a compact set K £ G Q 
such that k restricted to K £ is continuous while A(f2 \ lf e ) < e. Furthermore, there exists 
a continuous function k £ : R d x R d — >• R + with compact support such that k e = k on K e , 
supp A; e C fi e and ||fc e ||oo < C2. We estimate 

\\k-k £ \\ L i (n) = I \k-k £ \= I \k-k £ \+ I \k-k £ \<c 2 e. 
Jn J Ke Jn\K £ 

Define now k^(x, y) = C£([x] n , [y) n ) as above with k replaced by k £ . Then, since k £ is 
Riemann integrable with compact support, k e n converges in L^R 2 ^) for n — » oo to k £ . 
Moreover, for e < (4no)~ 1 , n > e~ x and by the definition of the conductivity functions 

\\K - k n \\ L i {n) < \\k £ - fe||£i(n e ) < \\k £ - k\\ L i {n) + \\k £ - fc||jr,i(n„-n) < c 2(l + ci)e 

Putting all this together we get for n large enough 

\\k - k n \\ L i {n) < \\k - k £ \\ L i {n) + \\k £ - k £ n \\ L i {n) + \\k £ n - fe„||n(n) < c 3 e. 

This directly yields (A5). Now Theorem 12. II applies. ■ 
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